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Abstract 

An effective integration method based on the classical solution to 
the Jacobi inversion problem, using Kleinian ultra-elliptic functions, 
is presented for quasi-periodic two-phase solutions of the focusing 
nonlinear Schrodinger equation. The two-phase solutions with real 
quasi-periods are known to form a two-dimensional torus, modulo a 
circle of complex phase factors, that can be expressed as a ratio of 
theta functions. In this paper, the two-phase solutions are explicitly 
parametrized in terms of the branch points on the genus-two Rie- 
mann surface of the theta functions. Simple formulas, in terms of the 
imaginary parts of the branch points, are obtained for the maximum 
modulus and the minimum modulus of the two-phase solution. 
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1 Introduction 

The focusing nonlinear Schrodinger (NLS) equation is 

ipt+Pxx + ^lpfp = 0, ( 1 ) 

where p{x,t) is a complex held exhibiting focusing, viz., modulationally un¬ 
stable, behavior. The NLS equation is a well-known and thoroughly studied 
soliton equation applicable to a wide variety of problems, in which there 
is a simple balance between dispersive and nonlinear effects, ranging from 
shallow water waves to optical communication systems, see PP and the ref¬ 
erences therein. The integration of multi-phase quasi-periodic solutions of 
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the NLS equation by algebro-geometric means has also been studied exten¬ 
sively [a El HI El E]. A crucial feature of the algebro-geometric technique is the 
determination of reality conditions which must be satished by the integrals 
of motion and the invariant spectral curve, in order for the N-phase quasi- 
periodic solution with real quasi-periods to satisfy the scalar NLS equation, 
not merely a complexihed coupled version of the NLS equation. Equiva¬ 
lently, the constants of integration in the algebro-geometric N-phase solution 
must satisfy reality conditions, so that the initial conditions of the solution 
are consistent with the integrals of motion used to construct the solution. 
The complete characterization of all reality conditions on the constants of 
integration in the multi-phase solution is a crucial feature in the correct de¬ 
termination of the solutions. Reality conditions for quasi-periodic solutions 
have been studied for the NLS equation [1] and for the sine-Gordon equa¬ 
tion mi, including very concrete formulas in the case of the sine-Gordon 
equation. 

The detailed study of two-phase solutions of the focusing NLS equa¬ 
tion Q is important for understanding the small-dispersion limit near the 
gradient catastrophe [9] and, in particular, the generation of rogue waves 
near the gradient catastrophe [10]. Special classes of ultra-elliptic solutions 
of vector NLS equations have also been of much interest recently mi [121 [131 
[T3|, including their modulation equations [T^- For elliptic solutions of the 
NLS equation, Kamchatnov miis] has demonstrated an effective integration 
method that clarihes the dependence of the solution on the branch points 
of the underlying elliptic curve. The purpose of this paper is to extend the 
effective integration method of elliptic solutions to two-phase quasi-periodic 
solutions, viz., the ultra-elliptic solutions. As is known HE], the real quasi- 
periodic solutions form a single two-dimensional torus in the Jacobi variety 
of the invariant spectral curve, modulo a circle of phase factors. In this 
paper, the two-phase solutions are integrated explicitly using classical meth¬ 
ods to solve the Jacobi inversion problem in terms of Kleinian ultra-elliptic 
functions mm and, consequently, in terms of Riemann theta functions ex¬ 
plicitly parametrized by the branch points of the genus-two Riemann surface. 
The resulting formulas are similar to previously known algebro-geometric for¬ 
mulas for the solutions [6], except that the dependence of the constants of 
integration on the branch points is more explicit. In particular, it is shown 
that, for all smooth two-phase solutions, there is a simple dependence of the 
maximum modulus and the minimum modulus of the two-phase solution on 
the imaginary parts of the branch points, see Theorems and 
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2 Lax Pair and Recursion Relations 

The integrability of the NLS equation ([^ is established through the equiva¬ 
lence of system ([^ and the commutation of a Lax pair of linear eigenvalue 
problems, 

V’x = W, = VV', (2) 

/ —iX ip \ 

y ip* i\ J ’ 

(3) 

/ —2iX^ + ijpj'^ 2iXp — px \ 
y 2iXp* + Px 2iA^ — 

p* denotes the complex conjugate of p, and A is the spectral parameter in 
the inverse spectral theory of the integrable system. 

The commutation of the Lax pair of linear operators (|^ implies the NLS 
equation in zero-curvature form 


where 

U = 

V = 


Ut-V, + [U,V] = 0, (4) 

where [A,B] = AB — BA is the usual operation of matrix commutation. The 
stationary N-phase solutions of the NLS equation are, by dehnition, solutions 
of the stationary NLS hierarchy defined by 

^x = [U,vI/], (5) 

in which the solution matrix T is polynomial in the spectral parameter A. 
The time-dependent N-phase solutions are then obtained by explicitly con¬ 
structing the compatible time-dependence of T such that 

T, = [V,T], (6) 


which, in turn, implies the zero-curvature representation of the NLS equa¬ 
tion Q. 

The commutation operators in equations (|^ and imply that the trace 
of T is constant and so, without loss of generality, constant multiples of the 
identity matrix may be added to T, and we may assume that the trace of 
T is zero. Furthermore, the commutation structure also implies that the 
characteristic polynomial of T is a constant, providing integrals of motion 
that enable the integration of the N-phase solutions. 

Therefore a solution to equation (|^ is constructed of the form. 




( ^11 ^12 \ 

V ^21 -^11 J ’ 


(7) 
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so that the Lax pair of equations ([^ and (@ become 


'Liix = 12+ ip'^ 21, 

^12^ = -2ip^ii-2iX^i2, 

'^21x = 2ip*^ii+2iX^21, ..X 

^iit = -{2iXp*+ pI)^i 2 +{2iXp-p^)^2i, 

'hi2t = -2{2iXp-p^)^ii + 2{-2iX‘^+ i\p\‘^)^i2, 

^2it = 2{2iXp*+pl)^ 11 +2{2iX^-i\p\^)^21- 

N-phase solutions correspond to 'L which are polynomial in A. Substitu¬ 
tion of the series 

OO 

^11 = 1+E/nA-", 

n=l 

OO 

'J'l2 = X] 9n^~^, (9) 

n=l 

OO 

flt21 = ^ hnX~'^, 

n=l 

into the stationary equation (|^, produces the Lenard recursion relations for 
the entries of d', 

9n = -Pfn-i, 

~ P* fri-l, (10) 

fn = ijphn- P*gn dx. 

It can be shown that the integrands are always exact 1201, SO that the entries 
of d' are differential polynomials in p and p*. The homogeneous solution, for 
which all the constants of integration are set equal to zero, is, up to 0(A“^), 


/l 

= 0, 

9i = 

-P, 

hi = 

-p*, 

/2 

= -¥p\", 

92 = 

l-dp 

2''9a:’ 

h2 = 

i-dp* 
2 ' dx 

/a 


93 = 

iS + ipIpt. 

hs = 

1 a^p* 
4 dx^ 


( 11 ) 


3 Two-phase solutions 

Two-phase solutions correspond to a third-degree polynomial solution for T, 
the most general such solution is 

Til = -iX^-ic 2 X^+ {\i\p\^-ici)X + \{pp*'-p'p*) + \ic 2 \p\'^-icii, 
Ti2 = ipX"^ + (-|p' + ic2p)X - lip” - lip\p\‘^ - \c2p' + i+p, 

^21 = ip*X‘^+ {\p*'+ ic2P*)X-\ip*”-\ip*\p\^+ \c2P*'+ icip\ 

( 12 ) 
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where Cq, Ci, C 2 G M are constants of integration. The stationary eqnation for 
this solntion is 


p"' = —6|p| V + 2ic2(p" + 2p|p|^) + Acip — Sicop. 


(13) 


The solntion T of eqnation (12) has real symmetry which mnst be satished in 
order to obtain two-phase solutions of the scalar NLS equation ([^, instead 
of solutions to a complexihed pair of coupled NLS equations. 


Theorem 1 (Reality Condition). 

^12(A) = -(T2i(A*))*, 


(14) 


hence the two roots pi and p 2 o/'h 2 i(A) = 0 are the complex conjugates of 
the two roots o/4/i2(A) = 0. In particular, 

Til = -iX^ - ic2)? + - ici)\ + \ii'2 + \ic2y1 - icQ, 

Ti 2 = ip{\- P2), (15) 

T 21 = ip*{\-p\){\-p*2), 

where z/i, z /2 G M are real variables, 


z/i = |p|2 > 0, 

1^2 = i{p*p' — pp*'). 


(16) 


The two roots pi,p 2 G C of the equation T 12 = 0 are analogous to the 
Dirichlet eigenvalues of the KdV spectral problem and, hence, are referred 
to as Dirichlet eigenvalues in this context. The solution p{x, t) of the NLS 
equation ([^ is recovered from the Dirichlet eigenvalues by the trace formulas. 

Lemma 1 (Trace Formulas). The Dirichlet eigenvalues pi and p 2 satisfy the 
trace formulas. 



Pi + P 2 — ~ C 2 , 

-1 // 1^-1 / 

(17) 


= -l^-lbP + l^C2|+Ci. 

Also 


1^2 — ~2z/i(/ii + P 2 + Pi + P 2 A- 2 C 2 ), 

(18) 

and 

p' = 2ip{pi + P2 + C 2 ), 

p" = ~‘2p{2pip2 vi 2 c2Pi -\- 2 C 2 P 2 + 2c 2 — 2ci). 

(19) 
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The evolution of the Dirichlet eigenvalues is governed by the Dubrovin 


equations obtained by substitution of equation ( [15| ) into the second and hfth 
equations of 


differentiating and substituting A = /xi or A = /i 2 . 

Lemma 2 (Dubrovin Equations). The Dirichlet eigenvalues /ii and fi 2 satisfy 
the system of Dubrovin equations, 


dfii 

0 ^ll(/4l) 

dx 

/41-/12 ’ 

9/42 

0 ^ll(/42) 

dx 

- ^ I 

M2-MI ’ 

9/41 

dt 

= -2{p2 + C2)^ 

dfL2 

dt 

= -2{pi + C2)^ 


( 20 ) 


The Dirichlet eigenvalues lie on trajectories in the complex plane deter¬ 


mined by the Dubrovin differential equations (20) and the initial conditions. 
However, not all initial conditions satisfy the reality condition. In other 
words, not all initial conditions are consistent with the integrals of motion 
of the Dubrovin equations and the assumption that the zeros of T 12 and 
'^21 are complex-conjugates of each other. The allowed initial conditions, 
viz., the allowed trajectories, are not known a priori (unlike the KdV case 
in which the Dirichlet eigenvalues must he on certain intervals of the real 
line determined by the spectral problem) and must be determined in order 
to construct two-phase solutions from the Dubrovin equations. 

The effective integration method of Kamchatnov [T[ [16] uses the fact 
that the Dubrovin equations show that the /i—trajectories of the Dirichlet 
eigenvalues form a real two-dimensional manifold parametrized by the two 
real variables vi and i> 2 - In fact, by Ending the algebraic dependence of the 
/i—trajectories on vi and 1 / 2 , the reality conditions can be satisfied explicitly. 
Moreover, the differential equations satisfied by vi and i >2 come from the 
flow and t—flow of the Lax pair for T by substitution of equations (15) 
into equations 


X 


Lemma 3. The real variables ui and V 2 satisfy the following system of equa¬ 
tions, 


dui 

dx 

du2 

dx 

dv\ 

dt 


= 2iz/i(/ri-l-/i2-/it-/ia), 

= 4iz/i(/it/it - T 1 T 2 ) - 2c2^, 


_ di'2 

dx '■ 


( 21 ) 


Sg- = 8!!/i((;i,- ri)|/i2p + (/12-/i;)|/iip)-4 c2|? 


Lemma 4. Ifp{x,t) is a smooth bounded quasi-periodic solution of the NTS 
equation 0, then vi{x,t) = \p{x,f)\‘^ must have relative extrema as a func¬ 
tion of X and t. Critical points o/z/i(a;, f) occur at distinguished values of the 
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Dirichlet eigenvalues. In particular, equations (21) imply that 


dui 

dx 


dui 

~di 


0 -v^ (i) = /ii, P -2 = lA, or (a) = fii. 


( 22 ) 


Notice that relative extrema of ui are the possible bounds for the oscilla¬ 
tions of a smooth quasi-periodic two-phase solution. These relative extrema 
can only occur at the distinguished values of the Dirichlet eigenvalues. 


Theorem 2. If p{x,t) is a smooth bounded two-phase solution of the NLS 
equation 0, then vi{x,t) = \p{x,t)\^ must oscillate on an interval of non- 
neqative values whose endpoints are relative extrema ofuii^x,!). If ui > 0 and 
hi 7 ^ h 2 , then oi has a relative extremum only z/pi = /i^ and /i 2 = pA- 


Proof. Explicit calculation using the Dubrovin equations 


and equa¬ 


tions ( 21 ) shows that if pi A p2 and (i) pi = and p 2 = lA (ii) hi = lA^ 
then 

However, if 7 ^ /i 2 and pi = P 2 and > 0, then 


dx'^ dt'^ dxdt 


-2Xv\\^u(h)? < 0 , 


(24) 


viz., there is a saddle point, instead of a relative extremum, at the critical 
point. □ 


It is worth noting that in order to obtain a similar result for iV— phase 
solutions with N > 3, the solution p of the NLS equation must be considered 
as a simultaneous solution of N evolutionary flows in the integrable NLS 
hierarchy of equations, because two conditions on the critical points of the 
X— and f—flows of the solution are insufficient to imply the condition of 
equation (22) on more than two Dirichlet eigenvalues. 


4 Invariant Characteristic Equation 


The characteristic equation of T is an invariant of both the x—flow and the 
t—flow of the solutions. As such, it provides the integrals of motion neces¬ 
sary to integrate the differential equations. In particular, for the two-phase 
solutions, we seek the explicit parametrization of the and p 2 trajectories. 
The characteristic equation of T, viz., the invariant spectral curve, as dehned 


by equations (15), is 


det{iwl — T) = —w'^ -f- ^(A) = 0, 


(25) 
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where 


*(A) = = fJ(A - A,), 


(26) 


i=l 


and it naturally defines a hyperelliptic Riemann surface of arithmetic genus 
two with points P = (A,tc(A)) G lying on a complex algebraic curve 
defined by equation ( [25| . The curve is a branched two-sheeted covering of 
the Riemann sphere with two points over the point at infinity; it is assumed 
that the curve is nonsingular, i.e., the branch points Aj, for i 
are distinct. Let the point cx)’*' be defined as the point over infinity with 
A = cx) and w{\) = A^ -I- O(A^), similarly cx“ is the point over A = cx with 
w{\) = —A^-I-O(A^). The curve admits the usual hyperelliptic involution 
corresponding to sheet interchange, 

i : Jp2 J^2i ^(A,tc(A)) = (A, —rc(A)) (27) 

and, because the coefficients of M{X) are real, an anti-holomorphic involution 

*(A,n;(A)) = (A*,n;(A)*), (28) 

where we use the same symbol * for the involution acting on points of 
as well as for complex conjugation of complex numbers. Note that the anti- 
holomorphic involution (28) leaves the sheets of the covering of the Riemann 
sphere unchanged, as can be seen by considering the action of * on points in 
the vicinity of viz., if A G M, then 

*(A,n;(A)) = (A*,n;(A)*) = (A,n;(A)). (29) 

The symmetry of the curve expressed in the existence of the anti-holomorphic 
involution * places reality conditions on the branch points of the integrals 
of motion of the x— and f—flows. 

Corollary 1 (Real Curve). must be a real algebraic curve, viz., the branch 
points are either real or come in complex-conjugate pairs. 


The Dubrovin equations (20) can now be integrated [2] on the hyperel¬ 
liptic Riemann surface. 


Lemma 5 (Dubrovin equations reprised). Eguation (26) implies 

Tii(pj) = (30) 

for j = 1,2. Therefore, the motion of the Dirichlet eigenvalues is defined on 
by the Dubrovin eguations {20), 


diJLi 

dx 

dx 

at 

dU2 

at 


_ 


2i 

2i 


; \/^(M2) 
U2-U1 ’ 


-2(/i2 -h C2)^, 
-2(/ii -F C2)^. 


( 31 ) 











Lemma 6. The x—flow and t—flow of fii and fi 2 on are linearized by the 
Abel mapping, 


/il ^2 

f- dfii j_ f dfi2 


mo V - f{o 
Ml M2 

r Ml<^Ml I f M2t^M2 

\/^(mi) 


MlO 


M20 


\/-^(M 2 ) 


4zt, 

2ia; — 4ic2t, 


(32) 


where /iio and /i 2 o constants of integration, viz., the initial values for /ii 
and fi 2 - 


It is important to remember that the values of /iio and fi 2 o are not arbi¬ 
trary but must satisfy the reality conditions of the integrals of motion and, 
hence, on the loci of pi and /i 2 by the characteristic equation (25). The 
correct determination of these initial conditions is one of the main goals of 
this paper. 


5 Integrals of motion and the Dirichlet eigen¬ 
values 

The symmetric polynomials of degree i, Sj, i = 1,..., 4, of pi, /i 2 , /i 2 , must 

satisfy the following reality condition because they must come in complex- 
conjugate pairs. 

Corollary 2 (Symmetric Polynomials of the Dirichlet Eigenvalues). 


and 


Si, S 2 , S 3 , S 4 G M. 

(33) 

S 4 = A 1 A 2 P ^ 0 . 

(34) 


Moreover, the algebraic constraints on the real loci of p,i and fi 2 can 
be made explicit. Using the invariant spectral polynomial (25) and dehni- 
tion (26), the symmetric polynomials of degree i, Aj, i = 1,...,6, of Aj, 
i = 1 ,... 6 , the branch points of the invariant spectral curve, can be written 
in terms of the symmetric polynomials of degree i, Sj, i = 1,... ,4, of the 
Dirichlet eigenvalues pi,/X 2 , Aj, A 2 - 
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Lemma 7 (Symmetric Polynomials of Branch Points). 

= — 2 C 2 , 

A 2 = 2 ci + C 2 , 

A 3 = — 2 ciC 2 — 2 co + ^ 1^2 + 2 c 2 i^i + 

A4 = -|c2l^2 - - CiZ/i + + 2 CoC 2 + cf + 1/1^2, 

As = -|z/iz/2 - ^C2l^i + Col^l + |ciI/2 - 2 coCi + CiC2l^l + I^lSs, 

Ae = \c2U1U2 - C2C0Z/1 - ^CoU2 + Icjuf + cl + Ai/2 


( 35 ) 


Eqnation (18) and the hrst three eqnations of (35) determine the three 


constants Cq, Ci, C2, in terms of the branch points. 
Lemma 8 (Constants of Integration). 

C 2 = -|Ai, 

Cl = ^A 2 - |A?, 

Co = —IA3 + IA2A1 — j^Af. 


(36) 


Eqnation (18) and the last three eqnations of (35) show that loci of the 


Dirichlet eigenvalnes can be parametrized in terms of and 1^2 by solving 
for Sj, i = 1,..., 4. 

Lemma 9 (Symmetric Polynomials of Dirichlet Eigenvalnes). 

Si = -it + Ai, 


E 2 

S 3 


= -|Aig - + 64 At+|AfA 2 IA1A3+A4 ^ ^ 1^2^ 

= i ^22 + (l|A 2 -lA 2 )^-|AiI/i 

, ffeAf-lA3A2 + iA2A3 + |AiA2-lA2A3+A5 , i , 
--^ 


(37) 


= |AiI/ 2 - J + i^l^2 - iA3)^ 


I'l 


+ + IA 1 A 3 . 

Note that Aj G M, i = 1,..., 6, are real parameters and z/i, z /2 G M are 
real variables. Thns the fonr roots /xi,/i 2 ,/r!,/i 2 G C of the qnartic eqnation 


JJ(/x - /Xi)(/i - n*) = jj!^ - Si/x^ + S 2 /i^ - S 3 /i + S 4 = 0, (38) 


2=1 


each lie on a two-real-dimensional manifold parametrized by z/i, 1/2 G M. The 
explicit solntions for pi and /X 2 can be fonnd nsing the standard procednre 
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for solving quartic equations. First the quartic polynomial in equation (38) 
is reduced by the transformation 


/i — p + -Si, 


so that equation (38) becomes 


+ S 2 /i^ 


Ss/i + S4 — 0, 


where 


5 2 

5 3 

5 4 


= -isf + s 2 , 

= |Sf — ^SiS2 + S3, 

= -JgSf+ jbS2S2-iSiS3 + S4 


(39) 


(40) 


(41) 


Note that Si G M, so the reality condition of Corollary [^transfers unchanged 
to the shifted variables S 2 , S 3 , S 4 G M, and S 4 > 0, because the roots of 


equation (40) still come in two complex-conjugate pairs. Completing the 


square of the hrst two terms of equation (40), 


(A" + tE2)" = i:3A+iS2-S4- 


If S 3 = 0, then equation (42) can be solved immediately to give. 


fi- -^\l -S 2 s 


Si 


4 S 4 , 


(42) 


(43) 


where s = ±1. The reality condition of Theorem implies that the shifted 
Dirichlet eigenvalues come in complex-conjugate pairs, so their values are 


constrained and the expressions for them in equation (43) can be simplihed. 


Corollary 3 (First Symmetric Polynomial Constraint). If S 3 = 0, then 
S 2 > —2'\/S4, and the explicit expressions for the Dirichlet eigenvalues fall 
into two subcases. 


1. //si < 4 S 4 , then 

hi = 
h 2 = 

2. //S 2 > 2 \/S 4 , then 


-S 2 -|- 2a^S 4 i h'ly ^2 + 2'\/S4, 


-S 2 -l- 2 '\/S 4 i ^iy f^2 T 2 '\/S 4 . 


(44) 


hi 


= ±^.^/S 2 + ^/Si-4S4, 


h 2 — 


±^n/S 2 


(45) 


Si 


4 S 4 . 
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If S 3 7 ^ 0, then a quantity z, to be determined, is added to the quantity 
in the parentheses on the left-hand side of equation (42), 

(/i^ -|- 2^2 4“ = 2 zfi ^ -\- ^^2 — ^4 + S22' -|- z '^. (46) 

The quantity z is chosen so that the right-hand side of the equation is a 
perfect square, which means that z should be chosen so that the discriminant 
of the right-hand side, a quadratic polynomial in fi, is zero, viz.. 


23 + E22" + (jE^-E4).--tEl = 0. 


(47) 


In fact the roots of equation (47) can be written down explicitly, if necessary. 


Lemma 10. //S3 7 ^ 0, then equation (4't) has at least one positive root. The 
three roots of the cubic equation ( 4 I) are given by the explicit formula, 

Q 


^ + YT^-Tv^ — - , -—, 

3 


(48) 


where 


(49) 


n —_Lv2 _ Iv 

V 36^^ 2 

p _ 1 ^3 Iv V _L 1 V2 

~ 216^2 + ig2-i3, 

A = Q^ + R\ 

and A <0 if and only if all three roots are real. A double root occurs if and 
only if A = 0. 

Definition 1. //S3 7 ^ 0, then z is chosen, without loss of generality, to be a 
positive root of Equation (4'^). 


Since z is a root of the cubic equation (|47|), the right-hand side of equa- 

(50) 


tion (46) is a perfect square, 

1 


Therefore 


(h^ + 2^2 + — 2z{fi + —S 3 )^. 


pf -|- s\/2zjl -|- 7-S2 -|- s — 5 ^S 3 -j- z — 0, 
2 


(51) 


where s = ±1. Finally, solving the quadratic equation (51) gives explicit 
expressions for the shifted Dirichlet eigenvalues, together with algebraic con¬ 
straints. 

Corollary 4 (Second Symmetric Polynomial Constraint). If E 3 7 ^ 0, then 
z > 0 is chosen to be a root of equation ( 4 I) and the reality condition on 
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the Dirichlet eigenvalues, viz., they must form two complex-conjugate pairs, 
implies that 

P- - ^2-^S3 + Z>0. 

The explicit expressions for fii = fvi — and fl 2 = 1^2 — in terms of ui 
and 1 Z 2 are 


(53) 


h - - \/f ± A 
P‘2 = \/f =*= + 2 :. 

In particular, iff^s 7 ^ 0, then fvi 7 ^ /i 2 - 

When S 3 7 ^ 0, 2 ; > 0, so the set of (i^i, z/ 2 ) for which S 3 7 ^ 0 and either 
p+ = 0 or p_ = 0 is equivalent to the set of solutions of 


2^(^ + S2)^ -S^ = 0. 


(54) 


Now 2 ; > 0 is also the root of the cubic equation (47). The set of points 
{^ 11 ^ 2 ) for which Equations (47) and (54) have a common root is given by 
the resultant of the two equations. The resultant of Equations (47) and (54) 
is a non-zero constant multiple of 


S2(16S4(S2 - 4 ^ 4 )^ - S^(4S2 - I 44 S 2 S 4 + 27S2)). 


(55) 


Theorem 3 (Boundary of the Real Parameter Set). Suppose ui > 0, so 
that S 2 ,S 3 , and S 4 are all defined, and also assume that S 3 7 ^ 0. Then 
the boundary of the set of points (z/i,z/ 2 ) satisfying the reality condition of 
Corollary^is a subset of the algebraic set 

16S4(S2 - 4S4)^ - S3(4S2 - I 44 S 2 S 4 + 27S^) = 0, (56) 


representing the set of (z/i,z/ 2 ) for which eguations (fl) and (Sf) possess 
a common root Zc- On this algebraic set, the common root Zc of the two 
eguations (fl) and (5f) is either 


or 


— + 2 V ^2 + I 2 S 4 > 0, 


(57) 





1 

3 



+ I2S4 < 0. 


(58) 


The effective integration of the pi and p 2 equations is accomplished by 
transforming to the evolution equations for ui and V 2 . The equations (21) for 
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Ui and 1^2 become, in the shifted variables (ii and /i 2 , 

^ + fi2 — fi*i — fi*2) ■, 

^ ~ A 1 A 2 ) + (-|Si - 2C2)^, 

dv\ dv 2 

dt dx ’ 


( 59 ) 


Sg- = 8 foi(((ii-ADImP + (A 2 -/‘ 2 )l/‘iP) 

+(-lEf - 2 c 2E, - 4ci)|a - (E, + 4c2)|f. 

Eqnivalently, the explicit expressions for fii and /i 2 can be nsed to write down 
the complicated bnt explicit eqnations for the evolntion of vi and z/ 2 . 


dvi 

dx 

du2 

dx 

dvi 

dt 

du2 

dt 


-2V2iy,{^+^), 

- ^/pZ) + y2ni(Si + 4 c2)(7^ + 
du2 

dx ’ 

-A^/2vl{{z + p-)^ +iz + p+)^) 
+V2i/i(-|S2 - 4 c2Si - 8c2)(ypT + ypT) 
-4z/i^(Si + 4 c2)(^^ - ^). 


Assnming that a two-dimensional region of allowed initial valnes for (i/i, 1 / 2 ) 
exists, it mnst be bonnded by a snbset of the algebraic set defined by eqna- 
tion (56). Moreover, assnming existence of solntions to eqnations (60), viz., 
a region in which both > 0 and p_ > 0 , the slopes of the bonndary cnrves 
can be determined directly, viz., if > 0 and p+ 7 ^ 0 , then p_ = 0 implies 


dv2 

dui 




1 

2 


El — 2 c2, 


(61) 


in which either x or t is held constant. Similarly, if > 0 and p_ 7 ^ 0, then 
p+ = 0 implies 

^ = Jsi- 2 c 2 . (62) 

oui 2 

Thns, for ui > 0, the one-dimensional varieties dehned by p+(i^i, 1 ^ 2 ) = 0 or 
p_(z/i,z/ 2 ) = 0 dehne smooth cnrves z /2 = z/ 2 (z/i) which cross transversely at 
points where 2 ; 7 ^ 0 and p+ = p_ = 0 . At points where z = 0 and p+ = p_ = 0 
the two cnrves p+ = 0 and p_ = 0 meet tangentially. 


6 Dependence of Extrema on the Branch Points 


Having obtained the explicit expressions for pi and p 2 in terms of ui and z /2 
in Corollary these expressions can be nsed to solve the z/i and z /2 eqna¬ 
tions of eqnation (21). Assnming that a smooth bonnded two-phase solntion 
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exists, there must be at least one positive relative extremum for i/i, which 
can occur only if /xi = /i^ and ^2 = /^25 provided /ii 7 ^ /r 2 - It is now possible 
to characterize explicitly all cases of the disposition of the branch points of 
J ^2 for which the solution for Vi could possibly oscillate on an interval of 
non-negative values. 


Lemma 11 . /xi = and ^2 = */ */ ^3 = 0 ^2 = — 2 \/S 4 . 

Also, fii = 1^2 or fii = /i 2 if and only z/ S 3 = 0 and S 2 = 2 \/S 4 . 


Proof. The explicit expressions in equations (44), (45), and 
stated results. 


imply the 

□ 


Lemma 12 . If p{x,t) is a smooth hounded two-phase solution of the NLS 
equation 0 and vi{x,t) = |p(a;, t)p, then the non-zero extrema of vi{x,t) 
must occur at points (z/i, 1Z2) where S3 = 0 and S2 ~ 4S4 = 0. 

Proof. Theorem 1^ shows that extrema can only occur when either (i) pi = p2 
or (ii) Pi = p\ and p2 = pi- Together with the previous lemma, the result 
follows. □ 


Lemma 13. Suppose the quantities vi,V 2 £ 1^ are replaced by complex vari¬ 
ables h'i,iz 2 E C in the rational functions S2, S3, and S4. Then S3 = 0 and 
S 2 ~ 4 S 4 = 0 if and only z/P(z/i) = 0, where Pivi) is a polynomial of degree 
ten in ui G C. 

Proof. P(z 2 i) = 0 upon elimination of 1^2 from the two equations S 3 = 0 and 

S2-4S4 = 0. □ 

The power of the effective integration method presented in this paper lies 
in the fact, stated in the following theorem, that the roots of the polynomial 
equation P(z^i) = 0 can be expressed explicitly in terms of the branch points 
of the Riemann surface 

Theorem 4. If p{x, t) is a smooth bounded two-phase solution of the NLS 
equation 0 , then ui{x,t) oscillates on an interval bounded either by zero or 
by positive roots of the polynomial equation P(i^i) = 0. The polynomial Pivi) 
is a resolvent polynomial for the factorization of the reduced form of M{X) 
into two cubic factors. In particular, the ten roots of P^ui) = 0, together with 
the corresponding z /2 values, are given explicitly in terms of the six branch 
points Aj, z = 1 ,..., 6 , of by the formulas 


Ui = 

-Ixl 

z/2 = 

-lxiX2, 
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where 


Xl — -^7r(l) + -^77(2) + -^7r(3) ~ -^7r(4) ~ -^7r(5) ~ -^7r(6); 

X2 = + Ajp, + A=,3, - - A=,„, 


(64) 


and ( 7 r(l), 7 r( 2 ), 7 r( 3 ), 7 r( 4 ), 7 r( 5 ), 7 r( 6 )) G Sq, the permutation group of order 
six. 


Proof. If yUi 7 ^ /i 2 , then extrema of Ui (or saddle points) occur when pi = 
p-i, P 2 = p 2 hi = pi) either case, although we are only con¬ 

cerned with extrema) Pii'i) = 0. For each of these ten roots, the sextic 
factors into the product of two cubics Qi and Q 2 , 


^(A) — —4/ii(A)^ d-z/i(A —/ii)^(A —/i2)^, 
= Qi(A)Q2(A), 


(65) 


where 

Ql = ^d'ii(A) -f iydd|((A — /ii)(A — /r 2 ); 

( 66 ) 

Q 2 = id'ii(A) - iyFr(A -/ii)(A -/i 2 ). 

Notice that, in case pi = fi 2 , it is still true that Pir'i) = 0, and there is still 
a similar factorization of ^ into cubic factors, so the stated formulas for ui 
and U 2 are still true. The zeros of each cubic factor form two non-overlapping 
sets, each comprised of three distinct branch points. Using Vieta’s formulas 
to express the coefficients of the cubic factors Qi and Q 2 in terms of the 
three roots of Qi = 0 and the three roots of Q 2 = 0 , explicitly solving for 
ui and i >2 from the coefficients of the quadratic and constant terms of Qi 
and Q 2 , and eliminating the terms depending on pi and /i 2 , leads to the 
stated formulas for z/i and 102 . Notice that both z/i and U 2 are invariant under 
the actions of re-ordering the hrst three or the last three summands, and 
also swapping the hrst three summands with the second three summands, so 
that the number of values of either ui or U 2 is 6!/(3!3!2!) = 10, as expected. 
In fact, explicit calculation shows that if ui = —a^, then the degree ten 
polynomial Pir'i) is the resolvent polynomial for the factorization into two 
cubic factors, one of the factors being A^ -|- a 2 A^ -|- OiA -|- Oq, of the reduced 
sextic polynomial of the sextic polynomial ^{X). The resolvent polynomial 
can be calculated by successive eliminations of variables from the equations 
obtained by equating coefficients of the expanded product of the desired 
factors equal to the polynomial to be factored | 2 I]. □ 

Theorem 5. If a smooth bounded two-phase solution of the NLS equation 0 
exists, then the values of ui and 1/2 at critical points of ui {x, t) are given by 


= -Ixl 

^2 = -lxiX2, 


(67) 
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where 


Xl — -^7r(l) + ^ttW) + -^7r(3) ~ -^7r(4) ~ -^7r(5) ~ 


W( 6 ) ? 


X 2 = + Ajp, + A=, 3 , - - A=,„, 


( 68 ) 


Moreover, if Hi ^2 M a critical point, then the following expressions of the 
second-order partial derivatives at the critical point have a similarly explicit 
dependence on the branch points, 

I'lxxl'ltt - I'lxt = 256l/i^ii(/ii)^ii(/i2) 

3 

n ('^7r(A:) '^7r(4)) ('^7r(/c) -^77(5)) ('^7r(/c) '^7r(6)): 

k=l 

(69) 


and 

T^lxx 


= %ivi 


aI'ii(/ai) —4'ii(;a2) 
/A 1 -/A 2 


2(A7 i-( 1) A7 i-( 4) ) ( ^77(1) ^77(5) ) ( A77(2) A77(4) ) ( A77(2) A77(5)) 

2(A77(4) Aj 7(4)) (A77(4) '^■7r(6)) ('^■7r(3) '^ 71 ( 4 )) ('^■7r(3) '^ 77 ( 6 )) 

2(A77(2) -^77(5)) (-^77(2) -^77(6)) ('^77(3) '^77(5)) (-^77(3) -^77(6)) 

2^i(A77(4) A 77(4) ) ( A77(2) '^77(5) ) (-^77(3) -^77(6))) 

where (7r(l), 7r(2), 7r(3), 7r(4), 7r(5), 7r(6)) G Sq, where Sq is the permutation 
group of order 6. 

Proof. The first part of the theorem is just a restatement of the result of the 
previous theorem, viz., by Lemmaat a critical point of ui{x,t) either (i) 
/ii = fil and fi 2 = p 2 (ii) Pi — P* 2 - either case Pivi) = 0 by Lemmas 10 
and 12. According to Theorem vi and V 2 are given by the stated formulas 
for one of the permutations of the indices. The expressions for the second- 
order partial derivatives at the critical point, assuming /xi ^ pL 2 , follow by 
direct substitution of (i) pi = /ij and /i 2 = P 2 or (ii) pi = /i^ into the 
expressions for the second-order partial derivatives, resulting in formulas in 
terms of Tii(/i). Then the explicit expressions for /xi and fi 2 in terms of ui 


and z /2 are substituted into \hii(/i), using equations (39) and (44). Finally 


the explicit expressions for ui and z /2 in terms of the branch points, given by 
equation (67), are used. □ 


7 Branch Point Reality Conditions 

Using Theorem]^ the reality condition of Corollary [T] on the algebraic curve 
can be strengthened. 

Lemma 14. If Xi, X 2 , A 3 , A 4 , A 5 , Xq G M, then Pi^i) = 0 has no positive roots. 
Proof. The result follows immediately from equation (63). □ 
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Lorillllcl 15 . If 1^2 ^ -^ 1 ) -^ 2 ) -^ 3 ) -^ 4 ) -^ 5 ) -^6 ^ 

then P(z 2 i) = 0 has no positive roots. 


are all distinct, and Ai = A 2 , 


Proof. The proof is by contradiction. If Pi^i) = 0 and z/i > 0, then equa¬ 


tion (63) implies, for example, that A 3 -|- A 4 = A 5 -|- Ae and, also, that 
A3A4 = AsAe, since 1^2 is real. However these two equalities imply that A3 
is equal to at least one of A 5 or Ae, contradicting the distinctness of the 
branch points, viz., the non-singularity assumption on □ 

Lemma 16 . If 1^2 G 1^, Ai, A2, A3, A4, A5, Ae G M are all distinct, and Ai = A2 
and A 3 = A 4 , then Pi^i) = 0 has no positive roots. 

Proof. The proof of this lemma is similar to the previous lemma. If > 0 


in equation (63) and the corresponding 102 G 


in equation (63), then either 

□ 


As = Ae or, otherwise, Ai equals one of A 5 or Ag. 

Theorem 6 (Reality Conditions for Bounded Two-Phase Solutions). If the 
NLS equation & has a smooth hounded two-phase solution, and JP 2 is non¬ 
singular, i.e., its six branch points are distinct, then the branch points form 
three complex-conjugate pairs. = 0 has exactly four distinct non¬ 

negative roots, these are the possible non-zero bounds of the two-phase solu¬ 


tion, and the corresponding values of 1^2 given by equation (63) are also real. 
In particular, if Xi = A 2 = ri-t-isi, A 3 = A 4 = r 2 -t-is 2 and As = Ag = rs-l-iss, 
with 0 < Si < S 2 < S 3 , then the four real roots o/P(z/i) = 0 are 


0 < ..p < ..p < ..p < p-', 


(71) 


where 


.p = 

u? = 

( 3 ) 

= 

( 4 ) 

iz) = 


(si -f- S 2 — S3)^, 
(Si — S 2 + S 3 )^, 
(~Si -p S 2 + S3)^, 
(Si -f- S 2 + S3)^, 


(72) 


and 


^2 — 

‘^2 — 

P' = 


4(si -P S2 - S3)(riSi -p r2S2 - r3S3), 
4(si - S2 -p S3)(riSi - r2S2 + r3S3), 
4(-si -P S 2 + S 3 )(-riSi -P r 2 S 2 -P r 3 S 3 ), 
4(si -P S2 + S3)(riSi -P r2S2 + r3S3). 


(73) 


Proof. The preceding lemmas eliminate any other possibilities for branch 


points satisfying the reality conditions. Theorem 4, equation (63), provides 
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the explicit expressions for the four real roots. The explicit expressions for 
z/i show that the four real values are distinct, otherwise at least one of si, S 2 , 
or S 3 must be zero, which contradicts the assumption that the branch points 
are distinct. □ 


Lemma 17. If a two-phase solution of the NLS equation 0 exists, then it 
is bounded, viz., there exists a positive number M G M such that i'i{x,tY + 
i' 2 {x, ty < M for all x and t. 


Proof. First we show that vi is bounded. 

1. If 1^1 —)■ 00 and z /2 = oivi) and Ai 7 ^ 0, then S 4 ~ < 0, which 

is impossible since S 4 > 0 . 


2. If z/i —)■ 00 and V 2 = o{i>i) and Ai = 0, then S 4 
impossible. 


Wui 


< 0 , also 


3. If 1^1 —)■ 00 and 1^2 = then S 2 ~ ~ “^Aiz/i + ^ 1^25 and 


z 0. Hence 


^ S 3 
S 2 i —-!=. + 2 ' 
^plz 


—- or ± 
4 


■|Aiz/i + 


(74) 


so that, in any case, at least one of pi or p 2 will become negative, which 
contradicts the reality condition of Corollary 

4. If 1/4 —)■ 00 and ^ 00 , then S 4 ~ < 0, which is also impossible. 


VI 


The preceding cases are exhaustive and show that vi is bounded. 
Consequently, if |z/ 2 | —)■ 00 , then 




which contradicts S 4 > 0. Hence \i' 2 \ is also bounded. 


□ 


Lemma 18. If a smooth two-phase solution of the NLS equation exists 
in a neighborhood of a point where i>i{x,t) = then ui = 0 is 

a relative maximum of h'i{x,t). Moreover, the values of the corresponding 
solutions G and (p 2 , o' 2 \/^{p 2 )) G of the Dubrovin 

equations on the Riemann surface JP 2 clI this point are real and are given by 
explicit formulas. 


hi 

h 2 


?’l(s2+S3)+»'2(si+S3)+^3(si+S2)-\/ -4(r,y) 

2(81+82+^3) 

r'l(s2+S3)+T'2(si+S3)+r’3(.5i+S2)+y/ A{r,^ 
2(81+82+83) 


( 75 ) 
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where 


A{f, ^ = 4siS2S3(si + S 2 + S 3 ) + sl{r2 - + s|(ri - 

+s|(ri - r 2 )^ + 2siS3(r2 - r3)(r2 - n) (76) 

+2s2S3(ri - r3)(ri - r2) + 2siS2(r3 - ri)(r3 - r2), 
and A{f, s) > 0. The hyperelliptic irrationalities are given by 

= (Ti\pii - XiWui - XsWui - A5I, ( 77 ) 


for i = 1,2, with the hyperelliptic sheets given by ai = —1 and a 2 = 1 . 
Proof. Direct substitution shows that when {ui, V 2 ) = (^i^\ 


2(si + S2 + 53)2 ^ 

where 

A{r, ^ = 4siS2S3(si + S2 + S3) + s?(r2 - r^)'^ + s|(ri - r^f 

+s|(ri - r2Y + 2siS3(r2 - r3)(r2 - ri) 

+2s2S3(ri - r3)(ri - r2) + 2siS2(r3 - ri)(r3 - ra) 

= 4siS2S3(si + S2 + S3) + (si(r2 - r3) + S3(ri - r2)f 
Psl{ri - r^y + 2s2S3(ri - r3)(ri - r2) 

+2siS2(r3 - ri)(r3 - r2) 

> 0 . 


(78) 


(79) 


The positivity of A{f, ^ follows from the fact that the first expression in 
equation (79) for ^4(7*,^ is symmetric with respect to interchange of the 
indices i = 1,2,3. Furthermore, if si,S 2 ,S 3 > 0 and ri < r 2 < r 3 , then the 
second expression for A is positive. Thus S 2 = —2^/^ and, hence, pi = /i^ 
and Pj 2 = P* 2 ^ with 


/ii 


r"! (S2+S3)+f’2 (si+S3i+fs (.51+S2)- 

2(si+S2+^3) 

r’l(s2+S3)+i'2(si+S3)+r3(si+S2)+.\/ A{r,^ 
2(si+S2+S3) 


(80) 


Obviously, /^i, /i 2 G and pi 7 ^ /i 2 , as expected. Moreover, the real symmetry 
of the hyperelliptic curve implies that 


^{Pi) — \Pi ~ Xi\^\p,i — X'i\^\p,i — Asl^, (81) 


since /Xj G M for z = 1,2. 

The explicit expressions of Theorem 5 can be used to obtain the result, 
using standard results of differential calculus. Equation (30) is also used, in 
order to determine the correct hyperelliptic sheets, ai = ±1 or <72 = ±1, for 
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dx^ dt^ 


_ _ 
dxdt) 


the points (pi, al^J and (p 2 , 0 ‘ 2 '\/'^(h 2 )) on the Riemann surface 
Equation ( [^ becomes 

256z/^d^ii(/ii)d'ii(p2) 

-256z/f (Ti(T 2 ^ M{^i) ^ ^(p 2 ), 

256(si + S 2 + S 3 )siS 2 S 3 |Ai 

0 . 


A^PlAi 


A 5 PIA 3 


-^5 I 


> 


(82) 


Thus (T 1 CT 2 < 0 and a 2 = —cti. 


Equation (70) becomes 




'(Kii(Mi) —'3i'll(/^2) 


Aii-M2 
si+H+H., 




— —32(si + S 2 + S 3 )siS 2 S 3 — 8 S 1 S 2 IA 1 — A 


* 12 


(83) 


- 8 S 1 S 3 IA 1 


\ * 12 
^ 5 ! 


8 S 2 S 3 IA 3 - A 


* 12 


< 0 . 


Moreover, at this relative maximum 

□ 


Hence there is a relative maximum 
(Ji = —1 and (72 = 1 . 

Lemma 19. If a smooth two-phase solution of the NLS equation 0 exists 
in a neighborhood of a point {x,t,i'i{x,t)) = (a;,t, for i = 2 or i = 3, 
then (x, t, is a saddle point of ui = z/i(x,t). 

Proof. Consider the case where i = 3, the case for i = 2 is similar. As in the 
previous lemma, 

Ai 


S 2 = 


2(—Si + S2 + 53)2 ’ 


(84) 


where Ai = A(r, —si, S 2 , S 3 ) and A(r, s) is given by equation (79). In this 
case the sign of Ai = A(r, —Si, S 2 , S 3 ) is indeterminate, and 


_ ?~l(^2+^3)+^2(-8l+g3)+?'3(-Sl+S2)-AR 

hi 2 (-Si+S2+S3) ’ 

_ ?'l(^2+Jj3)+^2(-8l+53)+?'3(-Sl+S2) + VAr 

^■2 2(-Si+S2+S3) 


(85) 


If Ai 7 ^ 0, then pi 7 ^ /i 2 and the formula of equation (69) is applicable. If 
Ai = 0 , then /ii = P 2 £ but all the expressions for the second-order 
partial derivatives of h'i{x,t) have hnite limits as p .2 —t pi with = Pi- So 
the formula of equation (69) is still valid in the limit. In particular, regardless 
of the value of Ai, 

256z/i2Tii(pi)Tii(/i2) 


dV _ (dPo\2 _ 

dx'^ df^ 1 dxdt i 


— —256(—Si -7 S 2 + S3)siS2S3|Ai — A 3 PIA 1 — A 5 PIA 3 — A 


* 12 


< 0, 


(86) 
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so there is a saddle point. 


□ 


Lemma 20. If a smooth two-phase solution of the NLS equation exists 
in a neighborhood of a point {x,t,iyi{x,t)) = then ui = i'i{t,x) 

has a saddle point at {t,x, when S 3 < Si + S 2 , but z/i(t, a;) has a relative 
minimum of ui = > 0 when S 3 > Si + S 2 . 


Proof. As before, consider the formula given by equation (69). When ui = 


u 


( 1 ) 


dx"^ 1 dxdt i 


256z/^d'ii(/ii)d'ii(/i2) 


= - 256 (si + S2 - S3 )siS2S3|Ai 




“ Asl' 
(87) 

If S 3 < Si + S 2 , then the above expression is negative, so there is a saddle 
point. If S 3 > S 1 + S 2 , then the above expression is positive, and equation (70) 
becomes 


dx^ 1 Ml-M 2 

= 32 (si + S2 — ■S3)siS2S3 — 8S1S2IA1 — A3P 
+8S1S3IA1 — A5P + 8S2S3IA3 — Asp 

= 8(si + S2)(S3 - Si)(s3 - S2)(S3 - Si - S2) - 8SiS2(ri - r2f 
+8s2S3(r2 - rs^ + 8 siS 3 (ri - r^f 

> 8(si + S2)(S3 - Si)(s3 - S2)(S3 - Si - S2) - 8SiS2(ri - r2)^ 
+8s2(si + S2)(r2 - rs f + 8si(si + S2)(ri - r^)"^ 

= 8 (Si + S 2 )(S 3 — Si)(s 3 — S 2 )(S 3 — Si — S 2 ) 

+ 8 (s 2 (r 2 - r 3 ) + si(ri - rs)^ 

> 0 , 


so there is a relative minimum. 


□ 


Lemma 21. If > 0, then the four real points ^ 2 * 0 ) ^ = 1) ■ ■ ■ )4, 
given in Theorem 6, satisfying S 3 = 0 and S 2 ~ 4 S 4 = 0, are singular points 
of the one-dimensional algebraic set defined by equation (56), 




16S4(S^ - 4S4)2 - - I 44 S 2 S 4 + 27S^) = 0. 


(89) 


In particular, at (ui, ^ 2 ) = \ ^ 2 ^), fori = 1, 2, and 3, z/S 2 = — 2 '\/S 4 < 0, 

then {i'i \ 1 ^ 2 ^) is a node, but i/S 2 = 2 \/S 4 > 0 , then \ is an isolated 
point. However, at ( 1 / 1 , 122 ) = z/ 2 "^^), S 2 = —2\/S4 < 0, and 

is always a node. 
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Proof. When z/i > 0, equation (89) defines an algebraic set given by a poly¬ 
nomial equation Q(z^i,r' 2 ) = 0. At (i^i,r' 2 ) = substitution and 

simplification, using a computer algebra system such as Maple, shows that 
Q 0 , 0 , Cf 1/2 0 , and 


Qui^iQU2i'2 QviV2 -^'5i'S 2'^3(^l A S 2 “ 1 “ S 3 ) Cl]^Cl2®3 ^ 2 ’ 

where A" > 0 is a positive constant and 

ai = (r2 - r^f + (s2 + ss)^ 

02 = (ri - rs)^ -1- (si + ss)^ 

% = (^2 - ri )2 + (si + S 2 )^. 


(90) 


(91) 


Clearly 01 , 02 , and 03 are all strictly positive. In Lemma p! 8 | it was shown 
that S 2 = — 2 S 4 < 0 at (z/i, 02 ) = Hence, at (i^i, U 2 ) = z^ 2 ^^) 


Q UlUlQ l/2l'2 


Q 


U1V2 


< 0 


(92) 


and there is always a node at this singular point. 

In the case of each of the other singular points (z/|*\ z/ 2 *^), i = 1,2,3, 


equations (90) and (91) remain the same, except that the sign in front of s* 


changes from positive to negative. Thus the sign of the expression on the 
right-hand side of equation (90) is always the opposite of the sign of S 2 - 


Therefore, when S 2 > 0, there is an isolated point, but when S 2 < 0, there 
is a node. □ 


Lemma 22. The common root Zr of equations [fl] and (Sf) that exists at 
each point of the algebraic set (56), provided S 3 7 ^ 0, has a positive limit 
z = —S 2 > 0 at the node {i'f \ where S 3 = 0 and S 2 = — 2 '\/S 4 < 0 . 

Moreover, the same point is a regular point of each of the two 

distinct algebraic sets p+ = 0 and p_ = 0 which cross transversely as subsets 


of the algebraic set (56) at the node at , z/. 


( 4 ) „( 4 )^ 
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shows that S2 < 0 at \ Since S3 = 0, also, at this 
point, equations (47) and (54) have two common roots, z = 0 and 2 : = —S 2 - 


Proof. Lemma 


Therefore, the common root Zc of equations (47) and (54) that exists at each 


point of the algebraic set (56), provided S 3 7 ^ 0, must approach either 0 or 


— S 2 . The two possible expressions for this common root, given in Theorem]^ 
approach —S 2 and |S 2 , respectively. Since S 2 < 0 at (z^i , 1 ^ 2 ^'^), the limiting 
value must be 2 ; = — S 2 > 0 . 

Direct calculation shows that, at (z/j^\ z/ 2 "^^). 


> 0 , 


2(<ii+s2+S3)4(('5l + '52 )(-Si + S3)(S2 + S3)(si S2 S3) 

+SiS2(ri - r2)^ -I- SiS3(ri - r^)"^ + S2S3(r2 - rs)"^) 


-‘ 3 , 1^2 2 (si+S 2 +S 3 ) 


(93) 
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and 


(94) 




^2 


y/z 


3 , 122 - 


So at least one of p±,v 2 is nonzero. Thns is a regnlar point of 

at least one of p+ = 0 or p_ =0. Since the same point is also a node of 
the algebraic set (56), in a neighborhood of which Zc > 0 and p+p- = 0 


on the algebraic set, the two transverse one-dimensional snbmanifolds of the 
algebraic set at the node mnst correspond to p+ = 0 and p_ = 0 . 

□ 


Theorem 7. The set of possible initial conditions for the Dirichlet eigenval¬ 
ues Pi and p 2 for smooth two-phase solutions p = p{x, t) of the NTS equa¬ 
tion 0 is parametrized by a single compact connected region o/{(z/i,z/ 2 ) G 
M^}. The maximum value of vi = \p{x,t)\^ on this region is (si -f S 2 + ss)^, 
and the minimum value of ui = |p(a;,f)p on this region is (si-f 52 — 83)^ > 0, 
if S 3 > Si + S 2 , or zero, if S 3 < si -f- S 2 . 


Proof. The previons lemma shows that p+ = 0 and p_ = 0 cross transversely 
at the point {i'i,V 2 ) = where z > 0. These two one-dimensional 

algebraic sets divide a neighborhood of (z/i, z/ 2 ) = {i'^\ into four regions, 
exactly one of which (to the left of {ui, U 2 ) = {r'i^\ 1 ^ 2 ^'^), since must be a 
relative maximum) contains an open set where p+ > 0, p_ > 0, and z > 0. 
Thus there exists a region of initial conditions for ui and z /2 and, hence, for 
the Dirichlet eigenvalues pi and p 2 , on which all the reality conditions are 
satisfied. Since the smooth bounded solution constructed in this region must 
reach a maximum value of z/i = \p{x,t)\‘^, and z/i = is the only such 
maximum possible, the closure of the connected region where p+ > 0 and 
p_ > 0 , to the left of (z/i,z/ 2 ) = (z/|^\z/ 2 ^^) and bounded by the algebraic 
set (56), must be the unique compact connected region consisting of all the 
permissible initial conditions of (z/i,z/ 2 ). This region parametrizes the per¬ 
missible initial conditions for pi and p 2 through the explicit expressions of 
the first and second symmetric polynomial constraints given in Corollaries 3 
and 4, respectively. 

Finally, all smooth two-phase solutions of the NTS equation (|^ are bounded, 
by Lemma [1^ and so they must achieve a minimum and a maximum value for 
z/i = \p{x,t)\‘^. Therefore all such solutions must oscillate between the maxi¬ 
mum value and the minimum value of z/i = \p{x,t)\‘^ implied by Theorem 6 
and Lemmas 17, 18, and 19. □ 


The simple dependence of the extrema of the modulus of the two-phase 
solution on the branch points is consistent with the analogous result for 
genus-zero and genus-one solutions [ 1 ] and numerical simulations of genus- 
two solutions [ini [19] . For higher-phase solutions, analogous results require 
the inclusion of higher-time flows in the integrable NLS hierarchy. The spatial 
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and temporal flows of the scalar NLS equation are not sufficient to span the 
higher-dimensional torus of three-phase or higher-phase solutions. In general, 
the sum of the imaginary parts of the branch points of the higher-genus 
Riemann surface will be an upper bound for the maximum modulus of the 
solution. 


8 Two-phase solutions of the NLS equation 

In this section, the smooth two-phase solutions are constructed for the NLS 
equation Q. The solutions for the Dirichlet eigenvalues of such solutions 
were shown, in the previous section, to be a single compact connected two- 
dimensional manifold. As is well-known, the Dubrovin equations for the 
Dirichlet eigenvalues linearize via the Abel map onto the Jacobi variety of 


Ait + di = 


2ix — Aic 2 t + d 2 = 


/^i 

r ^ _ 


d/J,2 


^2 


f ^ _ 

v'^dd) 




fi2dfl2 


where 


di = 


d2 = 


MiO 

I 

Ai 

Mio 

I 

Ai 


dfJ,! 


fJ-20 

+ / 

A 3 


^^2 


dfi2 


(96) 


\/Sdd) \d^dd)’ 




/^20 


(96) 




\/^dd) y^^dd) ’ 


where the constants of integration di and d 2 are not arbitrary but are deter¬ 
mined by an allowed pair of initial values pio and /i 2 o which satisfy, along 
with the initial values of the real parameters ui and z/ 2 , the algebraic con¬ 
straints imposed by the invariant algebraic curve. Moreover, this set of per¬ 
missible values for di and ^2 is partitioned into equivalence classes lying on 
distinct non-intersecting two-real dimensional planes in formed by trans¬ 


lation along the linear flows of equation (95) on the Jacobian variety of 


Solving equations (95) for pi and ^2 is the classical Jacobi inversion prob¬ 
lem. The details of solving the Jacobi inversion problem, in terms of the 
Kleinian elliptic functions a and ( defined on the genus-two Riemann surface 
are placed in the Appendix. If u' = ^ with 


u'l 


Ur, = 


SO that 


^^2 

r d/ii r d/12 

y \/^dd) \/*dd) 

A*! ^2 

r /iid/ii r /L2d/L2 

jd \d^dd) \/*dd) 


Ait + di, 

2ix 2iAit + (^ 2 , 


u' = iVX -I- iWt -I- d. 


(97) 


(98) 
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where 


W 

d 


(0,2f, 

(4,2A,f, 

(di, ^ 2 )'^, 


( 99 ) 


and di and ^2 are defined in equation (96), then equation (A-64) produces 


simple formulas for the symmetric polynomials of fxi and /X 2 - In particular, 


hi + /i 2 — 5 A 1 




du) 


d“'+/A° 


hlh2 


- “1^1 + ^^2+ ^i - ^In 


du) 


du) 




( 100 ) 


where (5i and 62 are constants of integration expressible in terms of functions 
of the branch points. The solution p to the NLS equation ([^ is obtained 
from the trace formulas in Equation (|T7|) and Equation (100[), viz.. 


|:lnp = 


2?(/ii + P 2 ) - *Ai 


In 

OU2 


dn'+f- 


du) 




du) 


— 2160 


fin 

OX 


d-^'+fA du) 


d-'+ir; 


du) 


- 2id2, 


atlnp = 


-4i/ii/i2 + 2iA2 - \iA\ + Aif Inp 


4if - In 

OU\ 


'^(“'+/aT 


I In 


d^'+f^ du) 


+ 2*Aif^ln 


Aidi — 2i52Ai. 


du) 

Au'+fA du) 


’'(“'+/a^ . 

Hence, the two-phase solution to the NLS equation ([^ is 

n(x t) = Jly ^(0 a{iVx+iWt+p+) 

P{X,L) a(d+) a(iVx+iWt+3-) 


Aidi — 2i52Ai 

( 101 ) 

( 102 ) 


where 


(103) 


= /aT +/aT +/aT 
r = + 

(j){x, t) = — 262X — ( 4 ( 5 i -|- 262Ai)t A 00 
(5i, (52 G C are given explicitly by equation (|A-66), and 0o £ is an arbitrary 


real constant. Note that (5i and 62 in 0 are, in general, complex because there 
is also a complex phase in the ratio of the a functions. These complex phases 
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cancel out, so that the solution remains bounded for all x and t. The values 
of z^i(0, 0) > 0 and (/xio, sio), (/i 2 o, S 20 ) G can be explicitly given in terms 
of the branch points of the curve by the formulas derived in the previous 
section. 

It is convenient to remove the exponential factors from the a functions in 
equation (102) and write the two-phase solution in terms of Riemann theta 
functions. The dehnitions of the theta functions and the constants in the 
following expressions are given in the Appendix. 

v(xt) = ( 104 ) 

V\X: I V UU.llle ’ 

where 


= (l){x,t) + {2aiiU + {ai2 + 0 . 21 ){2x+ 2Kit)) dui 

-\-(2a22(2x -\- 2Ait) (q(i 2 + Q)2i)4f) du2- 


(105) 


Lemma 23 . The exponential factor in equation (lOf) is a phase factor 

of modulus 1 , viz., 

f{x,t) = KiX + rL2t + (fo, (106) 

where the wavenumbers Ki,K2 G M are given by 


K, = 2(Ai + A^-iAi)-(a;-i)i2(^-^)-(a;-%(^-^), 

K2 = 4(-|Aip - |A^IA 2 ) - 2(ca“^)ii(^ - ^) - 2(a;“^)2i(^ - ^) 


0+ 

+Ai (2(Ai A^ - 4 Ai) - {uj~^)i2{% - ^) - (w“^)22(^ - ^)) , 


W+ 6 

.0+ e: 


and (fo eM. is an arbitrary phase. 


(107) 


Proof. Since is purely imaginary (see the Appendix for details), it is 
sufficient to show that the differences of logarithmic derivatives of 9 of the 
form, for j = 1,2, 

9: 


3 _ j 

9 + 9- 

are also purely imaginary. Consider the complex conjugate 


(108) 


du + J^^du) 


du + f,yu) 


= -yny du+ydu) 


■\uj du + ydu)+m + Tm', 


(109) 
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where the sum of the half-period characteristics is zero, i.e., an integer lattice 
translation for some m,m' G Similarly, 



^A 2 


r‘A 2 


du + 


du 


'A 3 


' 2 ^ 


-1 


du+ du\+m + Tm', 


'A3 


( 110 ) 

with the same m, m' G as previously, since we can integrate along the 
same paths in the calculation of the half-integer characteristics. 

Now using the fact that 6 is an even function and the partial derivatives 61 
and 62 are odd functions, and the transformation property of the logarithmic 
derivatives in equation (A-29), we obtain, for j = 1 , 2 , 


Hence 


e+\* 


e+ 

-) = 

—2nimt 

e+’ 

e-\* 

—27iim'j 

e~ 

^) = 

j 

e- ■ 






( 111 ) 


( 112 ) 

□ 


Finally, all the previous results can be summarized in the following theo¬ 
rem, in which the solution of the NFS equation ([^ is constructed, explicitly 
parametrized by the branch points of the Riemann surface with all re¬ 
ality conditions satished. Moreover simple formulas for the minimum and 
the maximum of the solution are determined, based solely on the imaginary 
parts of the branch points of 

Theorem 8 . The smooth two-phase solutions of the NTS equation 0 with 
non-singular invariant curve form a two-real-dimensional torus subman¬ 
ifold, modulo a circle of complex phase factors, of the Jacobi variety o/J^, 
with branch points Ai = Ag = ri - 1 -isi, A 3 = A 4 = r 2 -\-is 2 , A 5 = Ag = r 3 -\-is 3 , 
where r*, Sj G M, for i = 1, 2, 3, and Si, S 2 , S 3 > 0. The two-phase solution has 
maximum modulus Si -f S 2 + S 3 > 0. The minimum modulus of the solution 
is S 3 — Si — S 2 > 0, if S 3 > Si -|- S 2 , or zero, if S 3 < si -t- S 2 . In general, the 
solution has an explicit representation in terms of the 6 function associated 
with fls constructed in the Appendix, and the branch points of 1 (^ 2 , viz.. 


p{x,t) 

where 


(si + S 2 + 


0 ^(3 ) 0 {loo ^U~^{x,t)) 

e{lu-^(3+)e{lu-^u-{x,t))' 


U^{x,f) = iV{x — xq) -J iW{t — to) + 13^, 
U~{x,t) = iV{x — xo) + iW{t — to)-\-13~, 
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(113) 


(114) 




( 115 ) 


with V = (0, 2f, W = (4,4(ri + rj + rajf, 


where the initial values of the Dirichlet eigenvalues (/xi, and 

{h‘ 2-1 \/^ih' 2 )) given by eguations (80) and (81), the phase ^{x,t) is given 
by equation (10 (^ , and {xo,to) G is an arbitrary location for the maximum 
modulus of the solution. 


9 Example Solutions 

Two solutions are constructed using the formulas of Theorem]^ Both solu¬ 
tions are constructed using Maple and a simple integration algorithm, based 
on Simpson’s rule, to compute the necessary contour integrals on the Rie- 
mann surface For integrals involving points over inhnity or branch points. 
Maple’s algcurve package is used to compute a Puiseux expansion near the 
point of interest, which is integrated directly using Maple’s int command |23]. 
In Figure the branch points are chosen to be 

— 2 + i, —2 — i, —1 + i,—l — i,l + i,l — i, 

so that the imaginary parts of the branch points are Si = 1 ,S 2 = 

S 3 = 1. Thus the maximum of the solution’s modulus is si -f S 2 + S 3 = 

■S3 < Si + S2, so the minimum is 0 . 

In Figure]^ the branch points are chosen to be 

— 2 + 3i, —2 — 3i, —l + i,—l — i,l + i,l — i, 

so that the imaginary parts of the branch points are Si = 3, S 2 = 

S 3 = 1. Thus the maximum of the solution’s modulus is Si -t- S 2 + S 3 = 

S3 > Si -|- S2, so the minimum is S3 — Si — S2 = 1. 

10 Conclusion 

An explicit parametrization of all smooth two-phase solutions of the scalar 
cubic NFS equation Q is given in terms of the branch points of the invariant 
spectral curve and the loci of the Dirichlet eigenvalues. In particular, simple 
formulas are found, in terms of the imaginary parts of the branch points, for 
the maximum modulus and the minimum modulus of each smooth two-phase 
solution. Independently, a general theta function formula for the solution is 
given with parametrization in terms of the branch points of the genus-two 


(116) 

1 , and 
3. Also 

(117) 

1 , and 
5. Also 
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Figure 1: Modulus of a solution of the NFS equation Q with maximum 3 
and minimum 0. 



Figure 2: Modulus of a solution of the NFS equation Q with maximum 5 
and minimum 1. 
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Riemann surface and the values of the Dirichlet eigenvalues which satisfy the 
reality condtions. 

The dependence of the maximum modulus of the two-phase solution on 
the imaginary parts of the branch points is consistent with known results 
for zero-phase and one-phase solutions pQ and with numerical simulations of 
two-phase solutions [iDiiin]. For higher-phase solutions of the scalar cubic 
NFS equation ([^, the two-real dimensional linear flow on the Jacobi variety 
will not span the higher-dimensional real torus of solutions, as it does in the 
two-dimensional case, so the sum of the imaginary parts of the branch points 
will be an upper bound on the modulus of the quasi-periodic solution, not 
necessarily the maximum value. 

The parametrization obtained in this manuscript may be useful in charac¬ 
terizing modulations of two-phase solutions, for example, in the vicinity of a 
gradient catastrophe in which spikes are formed which are limits of two-phase 
solutions [91110]. More work is needed to extend the current results to the 
scalar defocusing cubic NFS equation and the Manakov system of coupled 
NFS equations [m US] 122], [23] . 

It is well-known (see and the references therein) how to obtain the 
solution of the NFS equation p{x, t) in terms of theta functions dehned on 
the Jacobi variety of the Riemann surface by examining the asymptotic 
expansion of the Baker-Akhiezer solution of the Fax pair. In that approach, 
constants of integration are chosen in a manner sufficient to satisfy the reality 
conditions. In this paper, necessary and sufficient conditions for the existence 
of smooth two-phase solutions are derived for all the integration constants, 
using the explicit representation of the /i-trajectories of the Dirichlet eigen¬ 
values, analogous to the results for elliptic solutions obtained by Kamachat- 
nov mus]. Consequently, a theta function representation of the solution 
p{x^ t) is also constructed, using the classical integration of the Abelian inte¬ 
grals by the Kleinian ultra-elliptic functions, with all constants of integration 
expressed in terms of the branch points of the Riemann surface. 


11 Appendix - The Jacobi Inversion Problem 


The inversion of the Abelian integrals in equation (32) for pi = pi{x,t) 
and p 2 = is called the Jacobi inversion problem. Classically [IT] the 

inversion problem for the symmetric polynomials pi + p 2 and piP 2 was solved 
in terms of the Kleinian sigma and zeta functions dehned on the odd-degree 
genus-two curves which generalize the more familiar elliptic functions, i.e., 
the Weierstrass sigma and zeta functions. A more modern treatment [T8] 
shows how the Kleinian elliptic functions can be used to solve the inversion 
problem for even-degree hyperelliptic curves such as J^ 2 - 
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11.1 Differential Identities on the Riemann Surface 


In order to solve the Jacobi inversion problem (32), several definitions and 
fundamental identities involving differentials on the Riemann surface are 
needed. The goal of the present manuscript is to keep technicalities to a min¬ 
imum and follow as closely as possible the classical approach of Baker mi, 
while using the more modern language found in [ 6 ] and [IE]. First a canonical 
basis of homology cycles, {oi, 02 ; bi, 62 }, is constructed on satisfying the 
intersection properties Oi o 02 = bi o 62 = 0 and Oi o bi = 02 o 62 = 1. The 
action of the natural hyperelliptic involution t on is extended in a natural 
way to the homology cycles, and the cycles are chosen such that 


i(ai) = -Oi, 

i[a2) = - 02 , 

i(bi) = -bi, 

i(b2) = -b2. 


(A-1) 


In particular, the arrangement of cycles is chosen as shown in Figure for 
the case of three complex-conjugate pairs of branch points which satisfy the 
reality conditions for the NTS equation ([^. Similarly, the natural action of 
the anti-holomorphic involution * on the cycles is 


*(ai) = -oi, 

*(02) = —02, 

*(bi) = bi-F 02 , 

*(b2) = b2 + 0i. 


(A- 2 ) 


Now consider the two holomorphic integrals that appear in the Jacobi 
inversion problem (32), 


dui = 


dX 
w ’ 


, XdX 
(iti2 — - 


w 


(A-3) 


The natural action of the anti-holomorphic involution on these two differen¬ 
tials is *dui = du*, for z = 1 , 2 , where du* denotes the complex conjugate 
of the differential. The periods of the above differentials around the basis 
cycles are defined, for i,j = 1 , 2 , as 



Notice that for the given choice of the basis cycles. 


(A-4) 



so that = 0 . 


(A-5) 
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Figure 3: Basis of homology cycles, genus two, for focusing NFS equation. 


A normalized basis of holomorphic differentials, canonically dual to the 
basis of homology cycles, can now be constructed. Let, for i = 1,2, 

2 

dvi = Y. \{^~^)ijduj, (A-6) 

i=i 

then the periods are normalized in the sense that, for i,j = 1, 2, 



(A-7) 


where 6 ij is the Kronecker delta symbol. A standard result, using Riemann’s 
bilinear relations, is that the determinant of the matrix u = (wp) is nonzero, 
and the matrix 

r = (A-8) 

is symmetric with a positive dehnite imaginary part. 

With the particular choice of the canonical cycles shown in Figure the 
symmetry of implies further information about the real part 3fJ(r) of r. 
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VIZ., 


= 


§dv* 

= fE 

bj k=l 

2 

= -E\{^~^)ik§ *dUk 

k=l bj 

2 

= - f duk 

k=l *(bj) 

2 

= - \{^~^)ik i’ duk 

k=\ 

= — E li^~^)ik{‘^^kj + 2a)fcj) 

k=l 

= —Tij — Sij, 


(A-9) 


where Oi = 02 and 02 = Oi, so that the hrst and second colnmns of u are the 
same as the second and first colnmns, respectively, of cn, and 


h=(? J)- (A-IO) 

Therefore 

»(r) = 4 ( ? J ) ■ 

We now introdnce two differentials of the second kind associated with dui 
and du 2 , see DU , viz.. 


^ 4A4-3AiA3 + 2A2A2-A3A^, 

CLTI — ---dA, 

4u? 


dro = 


2A3 - AiA^ 
Aw 


dX. (A-12) 


These differentials are holomorphic except for poles of the second kind at the 
points at inhnity, and they satisfy the following identity of two-differentials. 


d w + s 1 
dz X — z 2w 


dz dX + dui{X)dri{z) + du 2 {X)dr 2 {z) 


F{X, z) + 2ws dX dz 
4(A — z^ w s 


(A-13) 

where w‘^ = ^{X),s‘^ = ^{z), and the symmetric fnnction F(A, z) = F{z, A), 
where 


F{X,z) = 2Ae-A5{X+z)+Xz{2A4-A3{X+z))+X^z‘^{2A2-Ai{X+z))+2X^z^. 

(A-14) 


The periods of the differentials dri and dr 2 are defined, for i,j = 1,2, as 



(A-15) 
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The periods satisfy the following relations, 


/T 

OJU 

/ T 

— OJ U) = 



rjr] 

— rj rj = 

. .„/T 


UT] 

— U T] = 


(A-16) 


where denotes the transpose of the matrix w, and I is the 2x2 identity 
matrix. Equation (A-13) shows that the symmetric two-differential on the 
right-hand side of the equation gives Klein’s symmetric integral of the third 
kind, 


(A, a; z, h) = 



F(A, z) -\- 2ws dX dz 
4(A — zy w s ^ 


(A-17) 


with logarithmic inhnities of coefficients 1 and —1 respectively at A = x and 
X = b. Notice that the symmetry of the integrand implies that (A, a; z, b) = 
XXf {z, b; A, a). 


Definition 2. The normalized differential of the third kind having a simple 
pole of residue -|-1 at the location A = x, simple pole with residue —1 at the 
location X = b, and zero periods around the cycles Oi and 02 is denoted by 
dllz,b- The normalized differential of the second kind having a single pole of 
order two at X = z with coefficient 1 and zero periods around the cycles Oi 
and 02 is denoted by dT^- 


Notice that dT^ can be obtained from ^ by differentiation of the latter 
with respect to the parameter z. Thus an application of Riemann’s bilinear 
relations using dvi, i = 1,2, and dU^^b shows that 






dz 


(A-18) 


Another application of Riemann’s bilinear relations to dll\^a and dU^^b, demon¬ 
strates the following symmetry. 



(A-19) 


Since JfX (A, a; z, b) is an integral of the third kind, it can be written in terms 
of the normalized integral of the third kind and two linearly independent 
integrals of the first kind, leading to the following dehnition, in which the 
symmetry of and equation (A-19) are used. 

Definition 3. Let a be the 2x2 symmetric matrix defined by the identity 


XXL (A, a; z, b) = 


dU 


z,b 


2 2 

2 ^ ^ Olij 
i=l j=l 


dui 


duj. 


(A- 20 ) 
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Differentiation of the above identity with respect to local parameters A 
and ^ gives the following identity between two-differentials, 


2 2 

^(A, z) “h dX dz _ , ^, V—^ V—^ 1 / ■^ \ 1 / \ /A \ 

—-= dT^{X)dz - 2^^aijdui{\)duj{z). (A-21) 


4(A — zy w s 


i=l j=l 


If the roles of A and z are interchanged in equation (A-13), using the sym¬ 
metry of the two-form on the right-hand side, then equation (A-21), can be 
re-written as 


2 2 


(iMi(2;)dri(A)-|-(iM2(^)dr2(A)-|-2 EE aijdui{X)duj{z) = dTz{X)dz- 

i=l j=l 


d s + w dXdz 
dX z — X 2s 


(A-22) 

Integration in A of Equation (A-22) about the basis cycles and b^, for 
k = 1,2, gives two identities of the following differentials of z. 


2 2 

-2r]ikdui - 2rj2kdu2 + 4 X] Zl ^ik duj = 0, 

i=lj=l 
2 2 

-2r][j^dui - 2r]2f,du2 + 4 Z Z oiijw'ikduj = 27ri dvk 

i=lj=l 


(A-23) 


Since dui and du 2 are linearly independent, the hrst of the identities in 
Equation (A-23) implies 


—‘2r]ik + daiicuifc + 4a2in;2fc — 0, 

—‘2''f]2k + 4ai2Ci;ifc + 4a22i^2fc = 0 . 


(A-24) 


Since au = 021 , the preceding equations imply a key relation between the 
period matrices of the differentials {dui,du 2 } and {dri,dr 2 } which will be 
essential in dehning the Kleinian sigma function, viz.. 


1 

a = -Tju 


-1 


(A-25) 


11.2 Riemann Theta Functions 


The Kleinian sigma function will be dehned using the Riemann theta function 
associated with the period matrix r of the normalized differentials dvi and 
dV2. 

Definition 4. The theta function ofv G C^, associated with the period matrix 
T, is 


9{v) 


Z 


^27ri{vini^V2n2)-\-7ri{Tiin‘l^2Ti2nin2-\-T22rL2) 



D n2=—oo 
p 27 rmR+ 7 rmrn 


n=—oo 


(A-26) 
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with quasiperiodicity on the period lattice C^/(J 0t) given by 
9{v + m + rm') = 

where m, m' G 

Definition 5. The partial derivatives of 9 are denoted by 


(A-27) 


^i(^) = 

92{v) = ||(v). 


(A-28) 


Moreover, the logarithmic derivatives of 9 possess lattice transformations of 
the form, for j = 1,2, 


d\n9 , , d\n9, . 

-{v + m + rm ) = —2mmj H—^—(w). 


dvi 


dvi 


(A-29) 

Lemma 24. The choice of the canonical cycles obeying the conjugate re¬ 
lations (A-2) implies that the real part of t satisfies equation (A- 11) and, 
hence, it can be shown that 


Similarly, for i = 1,2, 


9{v)* = 9{v*). 

9i{v)* = 9i{v*), 
9^{v)* = 92 {v*). 


(A-30) 


(A-31) 


Theorem 9. The theta function is zero at odd half-integer periods, viz., if 
m, m' G T? and mm' = mim'i + m 2 m '2 is an odd integer, then 

9(-m H— rm') = 0. 

^2 2 ^ 

Proof. 

OO . 

^ ^7rimn-\-ni{n+^m')T{n-\-^rn')-^rn'Trn' 

n=—oo 

OO . 

n=—oo 

OO . 


n=—oo 

OO 






= e 


n=—00 
00 

E e 

n=—oo 
00 


— 2 nirnn+ 7 Timn-\- 7 ri{n-\-^rn')r(n-\-^m') — ^m'Tm' 


_ ^—i-Kvam' g7r2mn+7rz(n+^m^)T(n+|m^) — 


n=—00 

= e-'~'9{\m + \Tm'), 


(A-32) 
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where the summation index is renamed, hrstly, n —)■ —n and, secondly, n —)■ 
n — m', without change to the sum, which is over all possible integer pairs 
n e □ 


Definition 6 . The half-integer periods are denoted as 

1 
2 

where m, m' G 


m2 
mi m2 


m2 + T2imi 




Theorem 10. The fifteen integrals 

r 

^ (A-34) 

Jx, 

for i,j = 1 ,..., 6 , with i 7 ^ j, of the normalized holomorphic differentials 
dv = {dvi,dv 2 y on the dissected Riemann surface ^ constructed from the 
canonical homology cycles in Figure between the fifteen pairs of distinct 
branch points of the Riemann surface ^ are egual to fifteen distinct non-zero 
half-integer periods, as follows, 


Ai,A2 _ 

1 

0 0 



.yAljAs _ 

1 

-1 

1 


A1, A 4 — 

1 

-1 1 



2 

1 0 



2 

0 

-1 


2 

0 0 

AijAs _ 

1 

‘ -1 


0 


yX\,\Q _ 

1 

‘ -1 

0 ■ 


^A2,A3 _ 

1 

‘ -1 

1 


2 

0 


-1 

2 

1 

0 


2 

-1 - 

-1 

A2,A4 _ 

1 

■ -1 

1 



yX2,Xc, _ 

1 

■ -1 

0 


^A2,A6 _ 

1 

■ -1 0 ■ 


2 

-1 

0 


2 

-1 

-1 


2 

0 0 

A3,A4 _ 

1 

II 

0 

0 



yX3,X5 _ 

1 

‘ 0 

-1 ■ 


/^A3)A6 — 

1 

‘ 0 -1 



2 

0 1 


1 

2 

0 

0 

5 

2 

1 1 


A4,A5 _ 

1 

■ 0 - 

■ 



/^A4,Ag — 

1 

■ 0 

-1 ■ 



1 

0 

0 

1_ 



2 

0 - 

1 


5 

2 

1 

0 

5 

2 

1 1 



(A-35) 

Proof. By examining the dissection of the Riemann surface given by Fig¬ 
ure 1^ it is possible to integrate between any two branch points on the dis¬ 
sected Riemann surface while remaining on the lower sheet of the two-sheeted 
covering, without crossing any of the basis cycles. However the same integral 
can be performed by tracing the same path projected on the upper sheet for 
which the integrand is the same except for multiplication by —1 due to the 
action of the hyperelliptic involution on dv. By keeping track of the crossings 
of the homology cycles on the upper sheet (so as to remain on the dissected 
Riemann surface), the equality of the two integration procedures leads to the 
stated results. □ 
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Corollary 5. 

^ ^ e{v^^’^*) = ^ 0(v^3,A6) ^ = 0. 

(A-36) 

Lemma 25. If dv) is not identically zero, then it has precisely two 
simple zeros at \ = \i, A 3 . 


Note that d(0) certainly does not vanish identically for all d since 0(0) >0 
when ^ g ^ ^ cases 0 (f^^ dv) is not identically zero, since it 

is not zero at A = Ae- 

Lemma 26. If 6(j^^ dv) is not identically zero and Zi 7 ^ Z 2 , then d{J^^ dv — 
dv — dv) has precisely two simple zeros at X = zi, Z 2 G 

Lemma 27. If 6{J^^ dv) is not identically zero, then 6(z) = 0 if and only if 
A G such that 

+ / dv. 
j\3 

The set of z G C^/ (1 0 r) having this property is called the theta divisor 
0. The theta divisor is a one-complex-dimensional subvariety of the two- 
complex-dimensional period lattice / (I 0t). 

Lemma 28. Suppose 0(f^^ dv) is not identically zero and e ^ 0, then 

6 {J^^dv — e) is not identically zero and has precisely two zeros A = Zi,Z 2 . 
Moreover, up to addition of integer multiples of periods. 



e = 



(A-37) 


Thus, with the exception of the one-comp lex-dimensional variety e G 0, 
the points ( 2 : 1 , 81 ), ( 2 : 2 , 82 ) ^ that satisfy Equation ( A-37[ ) may be viewed 
as well-dehned functions of the independent variable e G C^/(J 0t). 


11.3 Kleinian elliptic functions 

Definition 7. The fundamental Kleinian sigma function a of u & is 

(j(u) = (A-38) 

with transformations on the period lattice C^/ (2u 0 2 u') given by, for r = 

1 , 2 , 

a(Ul + 2U}lr,U2 + 2U2r) = e2^lA«l+‘^lr) + 2,?2r(«2+t^2r)(j('.y) ^ 

a{ui + 2uj[^, U2 + 2u'^}j = e2^1.0i+<^ld+2»?h02+^h)a(M). ^ 
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In general, ifm = {mi,m2f and m' = {m[,m2r couples of integers 

and 

flm = 2um + 2u'm', Hm = 2rim + 2rj'm', (A-40) 


then 


<t(u + Ji „) = 


(A-41) 


Note that the sigma function can be multiplied by a constant factor inde¬ 
pendent of u without altering any of the results of this paper. The simplest 
normalization sufficient to accomplish the Jacobi inversion has been chosen. 
The definition of the Kleinian sigma function is motivated by the following 
identities between integrals on the Riemann surface. 

Firstly, Equation (A-20) implies 


(A, a; zi, bi) + (A, a; 2 ^ 2 , 62 ) - 








i=l j=l 

Secondly, the identity 


dui 


'bi 


duj -|- 


'62 


2:2,62 


duj 


(A-42) 


dl^zi,bi T 


dv - A"; dv - jll dv)e{jl dv - dv - r dv) 


^^22 ,62 


rh2 


log 


dv - /r; dv - dv)e{J^^ dv - dv - dv) 
follows from considering the function 
dilx, dv - dv - dv) 


obi 

IXi 


d(fxe dv - fx! dv - dv) 


62 

A3 


exp 


/ dYlzijbi 


dU 


22,62 


(A-43) 


which must be constant because it has no zeros and no poles and zero periods 
around all homology basis cycles. Therefore it is equal to its value when 
A = a. A similar conclusion follows even if only the periods around Oi and 02 
are zero and the periods around bi and 62 are constant (but not necessarily 
zero). 

The two identities in Equations (A-42) and (A-43) can be written in terms 
of the Kleinian sigma function (A-38), using the fact that a = as 


jr(A, a; zi, bi) -F a; Z 2 , ^ 2 ) = 

^(fx] du - du - du)a(f"^ du - du - du) (A-44) 

^ dv - du - du)a{£^ du - du - du) ’ 
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where du = {dui, du 2 )'^ and dv = {dvi, dv 2 )'^■ 

Now Equation (A-13), in which we interchange (A, a) and {z,h) using the 
symmetry of the expression, and Equation (A-17) imply 

Z2,b2) = 


fZi /‘A 


d s + w dX 


’ bi J a 


dX z — X 2s 


dz+ 


'•Z2 


'b2 J a 


d s + w dX 
dX z — X 2s 


dz + 


CZl 


''Z2 


dui + / dui 


'bi 


' i >2 


rzi 


f*Z2 


dU2 + 




'6i 


'62 


dri + 


dro 


(A-45) 


The identities in Equations (A-44) and (A-45) combine to give 

-A 


- Im ~ /a 3 duMlx, du - du - du) 


r '62 


log 


r*2:i 


^(fxe du ~ du - /a 3 du)a{J^^ du - du - du) 

r-A 


rb2 


d s + wdX , n 
dz+ / 


rzi 


’ bi J a 
rz2 


dX z — X 2s 

pX 


'62 


^^±!^^dz+ (A-46) 

dX z — X 2s 


dui -|- / dui 


'bi 


'b2 


dri + 


"Zl 


f22 


dU2 + 


du- 


'bi 


'62 


dro. 


Lemma 28 and the discussion following Equation (A-43) shows that in 


Equation (A-46), for u',u'' ^ 0, we can write 

PZl PZ2 


» 6 i 


f*b2 


U = 


du -f- 


du, u = 


du + 


du, 


' Ai 


'A 3 


'As 


consider zi and Z 2 as functions of u' and hi and 62 as functions of u”. Thus 
Equation (A-46) becomes, 


log 


cr 


(it du - u')(t{£ du - u") d s + w dX 


iXe 


f-a 

IXg 


r* A 

)Xq^'^ '-JXg 

d s + w dX 


cr(/^ du — u")a(f^ du — u') Jbi Ja d)X z X 2s 


dz-\- 


(A-47) 


'62 J a , 


dX z — X 2s 


dz -|- {u[ — u”) / dr I -|- (^2 — U 2 ) / dr 2 - 


Differentiation of Equation (A-47) with respect to u[ for i = 1, 2, produces 

rX 




du — u' ] -|- C 


' Ae 


, , , dzi d Si + w dX 

du — u = 


' Ae 


du', 


dX Zl — X 2si 


dZ2 

du', 


d S 2 + w dX 


(A-48) 


+ 


dri, 


dX Z 2 — A 2s2 

where the zeta functions Q and pij functions for i,j = 1 , 2 , are dehned by 

d d^ 
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Equation (A-41) shows that, for i,j = 1,2, and an integer translation 
across the period lattice of the Jacobi variety of 

Ci{u + rim) - Ci{u) = {Hm)i = 2rjiimi + 217 ^ 2^2 + + 2 ? 7 ' 2 m 2 (A-50) 

and 


Using 


we hnd that 


<pij(u -\- rim) r*p'(^)' 

_ r^i zdz I r^2 zdz 

“2 - Jai — + Ja3 


dzi _ 812:2 

aiti ^2—21' 

dzi _ 81 

du2 21 — 22 ' 


dZ 2 _ 8221 

9 mi 21—22- 

922 _ 82 

du 2 22—21 ■ 


(A-51) 

(A-52) 

(A-53) 


In Equation (A-48) we now make the change of replacing the points (zi, si), (z 2 , S 2 ) G 
= 4^2 with their corresponding points under the hyperelliptic involution, viz., 

(zi, —Si) and (z 2 , —S 2 ). Notice that u'^ is changed by this transformation to 

- u'^ + 2ujiimi + 2ui2m2 + + 2a;-2m2 (A-54) 

for some integers mi,m 2 ,m']^,m 2 . Making the corresponding change in the 
right-hand side of Equation (A-48) and using Equation (A-53), we obtain, 
for i = 1 , 2 , 


/ dvi + Ci 

J a 

where 


' \q 


du + u' ] --fi{X,Zi,Z 2 ) = Ci 


' \q 


du + u' ] - -fi{a,Zi,Z 2 ), 
(A-55) 


f ( \ ^ 2^ '1 — '»>(A-21-22) I 8i(2i-A-22) I 82 ( 22 -A- 2 i) 

Jl{A,Zi,Z2) - izi-\){zi-Z2) ^ {Z2-X){Z2-Z1)^ 

f2{X, Zl, Z2) — (^^^_X)(z2-\) (2l-A)(2l-22) (22-A)(22-2 i) • 


(A-56) 


By adding dvi + dvi to both sides of Equation (A-55), we obtain 


r'21 


f 22 


dvi + dri+ dvi + (i 


PZl 


'A 3 

/■22 


' Ae 


du + u' ] - -fi{X,Zi,Z 2 ) 


dvi + dri + Ci 


'A 3 


' Ae 


(A-57) 


du + u'] - -f,(a,Zi,Zu), 


in which the left-hand side is symmetric with respect to X,Zi,Z 2 , and the 
right-hand side is the value of the left-hand side when X = a. Consequently 
the left-hand side of Equation (A-57) is independent of X,zi, and Z 2 . There¬ 
fore, for z = 1 , 2 , 


0 


r 

'Ae 


f*Zl 


'*22 


du + u'] = Ci + Xf.(\zi,Z 2 ) - I dri- I dri- / dvi, (A-58) 


' Ai 


'A3 
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where Q is independent of A,Zi,Z2. Since a is an arbitrary point on /C2, it 
can be set equal to the branch point at a = Ag, giving 


r-zi 


'■22 


Ci ( I du + u'j — Ci +Zi, Z2) 


dr. 


dr. 


dr,, (A-59) 


' Ae 


' Ai 


'A 3 


Now Ci(0) = 0, being an odd function, so setting A = Ae, zi = Ai, and Z 2 = A 3 , 
shows that C, = 0, for i = 1 , 2 . 

Direct calculation shows that, as A —)■ 00 =*", the singularities in the terms 
on the right-hand side of equation (A-59) cancel out, so that 

1 /-A 


A— >Ooi 2 


' Ae 


where 


11 — -1-2V ^1^2 8 2 

I2 = ±\{Zi + Z2 - \^l) + , 


= it, 

(A-60) 


(A-61) 


and and 82 are constants independent of Zi and Z 2 . Thus equation (A-59) 
becomes 


0 


du + u' \ =7, 


± 




'■22 


dr. 


dr,. 


(A-62) 


' Ag 


'Ai 


'A 3 


Equation (A-62), implies, for i = 1,2, 


Ci{u'+ du]-Ci{u'+ dw = 7+ - 7, 


(A-63) 


' Ae 


' Ae 


Substituting the explicit expressions for 'yf, we obtain the solution to the 
Jacobi inversion problem. 


Cl + /aT “ Cl (u' + /a 7 duj = -Z1Z2 - I A? -F \k2 + ^1, 
C2 (u' + du) - C2 (u' + /“ dv) = Z1A-Z2- |Ai -F 82, 


(A-64) 

where the constants 81 and 82 are independent of Zi and Z 2 , and so can be 
obtained from equation (A-64) by setting zi = (Ai,0) and Z 2 = (A 2 , 0 ) = 
(A^,0). Recall aij for i,j = 1 , 2 , is given by 


a = -rju ^ = 


an ai2 
021 0:22 


(A-65) 
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Hence 


(/a 6 + I\3 (fxe 


So = <5^ 


where 


— -(-|Aip - |A^ + IA2) + 2 an f^- dui + (0:12 + 0:21) J^- 


2^'^ yiiVe+ g 

S2 = ~(Ai + A^ — |Ai) 


+*^2 (^/a 7 du + - C2 (^/a 7 du + du 

= ~(Ai + A^ — ^Ai) + (q;i2 + 0 : 21 ) f^- dui + 2q;22 f^- du2 

+^(a;-^)i2(S-^) + i(a;-%(^-^), 

(A- 66 ) 


0 + = 
0 - = 

= 

= 

Ot = 

= 


A 2 

A 3 


A 2 

A 3 


H/aT^C^W + /; 

o{i--HCdu+j;:: 

du+j. 

e, (kH/aT' 

^2 dU + J; 

^2 (^c^"H/aT ^“ + /; 


'^3 

A 2 

A 3 

A 2 

A 3 


(A-67) 


Explicit expressions for zi and Z 2 can then be found from the quadratic 
formula for the roots of the quadratic equation 


z‘^ — (zi + Z2)z + Z1Z2 — 0. 


(A- 68 ) 


The corresponding values of si and S 2 for {zi, si), {z 2 , S 2 ) G ^ are, for 
A; = 1,2, 


Sk = i P22 i u + / du 
JXe 

PO 

+P12 {u'+ du 

Jxe 


P22\U + j du\ \ Zk 
Jxa 

Pi2 du 

Jxe 


(A-69) 
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